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Jet substructure at HERA* 

Claudia Glasman a (On behalf of the HI and ZEUS Collaborations) 

a Universidad Autonoma de Madrid, 

Departamento de Ffsica Teorica, Facultad de Ciencias 

Cantoblanco, E-28049 Madrid, Spain 

A review is presented of jet substructure measurements from the HI and ZEUS Collaborations at HERA. The 
results presented include tests of perturbative QCD, comparison of the properties of quark and gluon jets, the 
comparison of the pattern of QCD radiation in different hard scattering processes, determination of the strong 
coupling and studies of the underlying subprocesses dynamics and of the pattern of parton radiation. 



1. Introduction 

Perturbative calculations deal only with par- 
tons in the final state. However, hadrons and 
not partons are observed in the detector. The 
hadrons are the result of the fragmentation of the 
partons such that all the hadrons originating from 
a given parton are contained in a narrow region 
around the direction of the original parton, what 
is commonly called a jet. Therefore, the first step 
in comparing data and theory is to reconstruct 
the parton topology. This is best done with jet 
algorithms, which define the jets in a very spe- 
cific way. The jets obtained display the following 
characteristics: (i) the longitudinal momentum 
(Pl) distribution of the hadrons in a jet should 
scale with the jet energy; (ii) the transverse mo- 
mentum (pt) distribution of the hadrons should 
have a mean value of about 300 MeV. Thus, the 
mean angle (~ Pt/pl) between a hadron and the 
jet axis and the size of a cone which contains a 
constant fraction of energy will decrease with in- 
creasing jet energy. At high energies, this picture 
is modified since gluon emission from the original 
parton overcomes the fragmentation effects and 
so the structure of the jets is driven by radiation, 
which has the advantage of being calculable in 
perturbative QCD (pQCD). 

Thus, the investigation of the internal structure 



of jets gives insight into the transition between 
the parton produced in the hard process and the 
experimentally observed jet of hadrons. QCD 
predicts that (i) the jet substructure is driven by 
gluon emission off primary partons at sufficiently 
high jet transverse energy, E^ , where fragmenta- 
tion effects become negligible, (ii) jets originating 
from gluons are broader than those originating 
from quarks since the gluon-initiated jets radiate 
more due to their larger colour charge and (Hi) 
the jet substructure depends mainly on the type 
of primary parton (quark or gluon) from which 
it originated and to a lesser extent on the hard 
scattering process. The experimental results from 
HERA which confirm these predictions are pre- 
sented in this report. 

The jet substructure can be studied in terms of 
two different observables, namely the jet shape [T] 
and the subjet multiplicity |2J. The jet shape is 
defined as the fraction of the jet transverse energy 
that lies inside a cone in the n ~ 4> plane of radius 
r, Ex(r), concentric with the jet axis. The mean 
integrated jet shape is given by 
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To compute the jet shape, the jets can be recon- 
structed either with the cone [3] or the kx [I] algo- 
rithms. In contrast, the subjet multiplicity can be 
computed using only the jets reconstructed with 
the ItT algorithm and is obtained by reapplying 
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the algorithm until for every pair of particles in- 
side the jet, the quantity 

dij = mm(E Tl ,E Tj ) 2 [(i]i - n^) 2 + (<fo - fij) 2 } 

is above y cu t • {E^ ) 2 , where y cut is the resolution 
parameter. The mean subjet multiplicity is given 
by 

<n s bj(j/cut)) = JT— ^ "sbj(ycut)- 
jets i=1 

The calculations of the jet substructure us- 
ing QCD-based leading-order (LO) Monte Carlo 
(MC) models, such as Pythia [5], Herwig [5J, 
Ariadne [7] and Lepto [5], are approximations 
based on a parton-shower approach. In fixed- 
order QCD calculations, the lowest order pro- 
vides no structure for the jets since there is only 
one parton inside a jet. Higher-order terms give 
the lowest non-trivial contribution to the jet sub- 
structure. In particular, it is possible to obtain 
next-to-leading-order (NLO) calculations for jet 
substructure in neutral current (NC) deep inelas- 
tic scattering (DIS) in the laboratory frame from 
0(a 2 ) predictions since, in this case, three par- 
tons can be inside one jet. Thus, measurements of 
jet substructure provide a stringent test of pQCD 
directly beyond LO. The pQCD calculations of 
the jet shape and subjet multiplicity are obtained 
via 

(1-VW> = 

j" dE T (E T /E^)[da(ep -> 2partons)/dE T ] 

Ojet 

Kbj(ycut)) = 
1 x 

H ( j - 1 ) • (7 s b j , j (t/cut ) = l+Ci a s +C2a 2 s . 

2. Tests of pQCD using jet substructure 

At HERA, the main source of jets is photopro- 
duction. At LO QCD, two processes contribute 
to the jet production cross section, resolved, in 
which the photon interacts via its partonic struc- 
ture, and direct processes, in which the photon 



interacts as a point-like particle. Both processes 
give rise to two jets in the final state. The observ- 
able x° bs = (l/S 7 )(X; je t s ^T te " I ' j0t ) measures 
the energy invested by the photon in producing 
the dijet system and can be used to separate re- 
solved and direct processes since they populate 
different regions of phase space. Resolved pro- 
cesses give rise to jets of quarks and gluons in the 
final state. On the other hand, direct processes 
are dominated by quark jets, so the 7y jet depen- 
dence of the jet substructure is expected to show 
quark-like jets in the rear direction and gluon- 
like jets in the proton direction due to HERA 
dynamics, since for the dominant resolved dia- 
gram (g p q^ — > gq), the gluon goes forward. On 
the other hand, the jets in an inclusive-jet sam- 
ple of NC or charged current (CC) DIS processes 
are dominated by quark jets so that no signifi- 
cant dependence with rf ct is expected. Only for 
a dijet sample, which contains a larger fraction 
of gluon jets, some dependence with if et could be 
observed. 

The mean integrated jet shape was mea- 
sured [9 for inclusive-jet photoproduction us- 
ing the iterative cone algorithm with cone radius 
R = 1. Figure [1] shows the mean integrated jet 
shape as a function of r for different regions of 
rf ct . The jets become broader as if ct increases. 
The comparison with the QCD-based MC calcu- 
lations shows that models with only fragmenta- 
tion predict jets that are too narrow compared 
to the data. In contrast, models which include 
initial- and final-state QCD radiation describe the 
data well for — I < rf ct < 1. Therefore, parton 
radiation is the dominant mechanism responsible 
for the jet shape. 

The comparison of the measurements with pre- 
dictions for samples of pure quark- or gluon- 
initiated jets (see Fig. [2]) shows that the measured 
jets are quark- like for — 1 < rf ct < and become 
increasingly gluon-like as 77 JGt increases, as pre- 
dicted by QCD. 

To see more clearly the dependence of the jet 
shape with rj ]ct and E^ , the mean integrated jet 
shape as a function of if Qt and Ej? for a fixed 
value of r = 0.5 [10] is presented in Fig. [3] This 
analysis was done using the kx algorithm to re- 
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Figure 1. Mean integrated jet shape in inclusive- 
jet photoproduction. 
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Figure 2. Mean integrated jet shape in inclusive- 
jet photoproduction. 



construct the jets. The measured jet shape de- 
creases as rf et increases, which indicates that the 
jets become broader, and the measured jet shape 
increases as E 3 ^ increases, which indicates that 
the jets become narrower as E 3 ^ increases. The 
comparison with the predictions for quark and 
gluon jets shows that the broadening of the jets as 
rf et increases is consistent with an increasing frac- 
tion of gluon jets. The comparison with the pre- 
dictions for resolved and direct processes shows 
that the events are dominated by resolved pro- 
cesses at low Ej, , whereas the direct contribution 
becomes increasingly more important as Ej, in- 
creases. 

This last result was studied in more detail [11] 
by separating a dijet photoproduction sample 
in x° hs < 0.75, dominated by resolved pro- 
cesses, and x° bs > 0.75, dominated by direct 
processes. The jets in the resolved-enriched sam- 
ple are broader than those for the direct-enriched 
sample (see Fig. [4]). The QCD-based MC predic- 
tions give a good description of the data and show 
that the data is consistent with being dominated 
by resolved (direct) processes for i° bs < 0.75 
(if > 0.75). 

The mean integrated jet shape was also mea- 
sured in NC DIS in the Breit frame for different 
regions of jet pseudorapidity, rf^, and transverse 



energy, -E^? B , for a sample of dijet events [T2"] . 
The measurements show that the jet shape for a 
fixed value of r = 0.5 decreases as a function of 
77g*. The effect is more pronounced at low E^ B 
(see Fig. [5]). The data are well described by the 
predictions of the QCD-based MC models; this 
comparison shows that the observed jet substruc- 
ture is compatible with that of quark-initiated 
jets. 

Figures [5] and [7] show the integrated jet shape 
measured in NC DIS in the laboratory frame for 
a sample of inclusive jets as a function of r in 
different regions of ^ ct and Ep [TU]. The mea- 
surements are compared to NLO calculations [13] , 
which give a good description of the data for r > 
0.1. The if et and E^ dependence of the mean in- 
tegrated jet shape was measured for a fixed value 
of r = 0.5 (see Fig. [5J. No significant variation 
with rf ct is observed in this case, consistent with 
the sample being dominated by quark jets in all 
the phase space measured. It is observed that the 
jets become narrower as E 3 ^ increases. The NLO 
calculations give a good description of the data 
and show a sensitivity to the value of ot s . From 
the measured jet shape for E^ > 21 GeV, a value 
of a s (Mz) was extracted, a s (Mz) = 0.1176 ± 
0.0009 (stat.) ™e (exp.) (th.), which is 

consistent with other determinations from more 
inclusive channels and has a small experimen- 
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Figure 3. Mean integrated jet shape in inclusive- 
jet photoproduction. 



tal uncertainty. The theoretical uncertainties are 
dominant, especially that arising from higher or- 
ders; a better determination of a s from this ob- 
servable would be possible once higher-order cor- 
rections for NC DIS become available. 

The mean subjet multiplicity has been mea- 
sured in inclusive-jet photoproduction [TO , di- 
jets in NC DIS [12] and inclusive-jets in NC 
DIS [H] (see Fig. [5]). The conclusions are sim- 
ilar to those from the integrated jet shape. In 
addition, a value of a s (Mz) was extracted from 
this observable in NC DIS: a s (M z ) = 0.1187 ± 
0.0017 (stat.) ±° : °°j£ (exp.) t°;°°°g (th.) 



dijet sample 
A [ x° b « < 0.75 



C. Glasman 
dijet sample 




S 1 


x° b ' > 0.75 


> 

V 




0.8 




0.6 




0.4 


j J • H1 Data (prel.) 


0.2 


I j — Pythla 

" 7 direct 

resolved 





i i i i i i i i i i i i i i i i i i i 



0.2 0.4 0.6 0.8 1 

r/R 



0.2 0.4 0.6 0.8 



r/R 



dijet sample 




Figure 4. Mean integrated jet shape in dijet pho- 
toproduction. 
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Figure 5. Mean integrated jet shape in dijet NC 
DIS in the Brcit frame. 



3. Jet properties in different hard scatter- 
ing processes 

The measurements of the integrated jet shape 
in photoproduction and NC DIS were compared 
in order to study the properties of jets in different 
hard scattering processes [15] . Figure [TUk shows 
the comparison of the jet shape in inclusive- 
jet NC DIS and dijets in photoproduction with 
x° bs ^ 0.75 as a function of r. It is observed 
that the jets in NC DIS, which are dominated 
by quark jets, are narrower than those in re- 
solved photoproduction, which include a larger 
fraction of gluon jets, and similar to those in di- 
rect photoproduction, which are also dominated 
by quark jets. Figure [TOb shows the integrated 



jet shape for r — 0.5 as a function of 77 Jct for sam- 
ples of inclusive jets in photoproduction, dijets 
with i° bs > 0.75 in photoproduction and inclu- 
sive jets in NC DIS. The inclusive-jet photopro- 
duction sample shows a behaviour which is con- 
sistent with an increasing fraction of gluon jets as 
if et increases, whereas the dijets with i° bs > 0.75 
and the NC DIS samples, which come predomi- 
nantly from quark jets, show no significant de- 
pendence with rf et . 

The parton composition of the final state was 
studied in more detail by comparing the ?y jet and 
Ej? dependence of samples of inclusive jets in 
photoproduction and NC DIS with the predic- 
tions of pure quarks and gluons [TO] (see Fig. [TT]) . 
The NC DIS sample is consistent with being dom- 
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Figure 7. Mean integrated jet shape in inclusive- 
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inated by quark jets in the whole phase-space re- 
gion. The jets in the photoproduction sample are 
similar to those of NC DIS for low 77-i ct and be- 
come increasingly broader as ?? Jct increases, which 
is consistent with the increase of the fraction of 
gluon jets. For the dependence, the NC DIS 
jets are well described by the prediction of quark 
jets in the whole measured range. The jets in the 
photoproduction sample get narrower somewhat 
faster than those in NC DIS. 

The jets in the NC DIS sample were also com- 
pared to those of an inclusive-jet sample in CC 
DIS J6], also dominated by quark jets, so a sim- 
ilar dependence with e.g. E^ is expected. Fig- 



■ ZEUS 98-00 DIS 
— DISENT 



E J T ct >17GeV 
Hadron level 




20 30 40 50 60 70 

E*' (GeV) 



Figure 8. Mean integrated jet shape in inclusive- 
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ure [12] shows the mean subjet multiplicity for a 
fixed value of y cu t as functions of Ej! * and Q 2 in 
NC and CC DIS. The jets in both samples get 
narrower as E^ increases, but the jets in CC 
seem to be slightly narrower than those in NC, in 
agreement with the NLO predictions [T7]. These 
differences can be understood in terms of the dif- 
ferent Q 2 spectra in NC and CC DIS. 

The measurements in inclusive jets in NC 
DIS |15] were compared to similar measurements 
from CDF [18], D0 [H| and OPAL [20]. The jets 
in NC DIS and e + e~ are similar (see Fig. IT3]) . 
which is understood in terms of a similar par- 
ton composition of the final state, dominated by 
quark jets. The jets in NC DIS and e + e~ are 
narrower than those in pp collisions, which can 
be understood from the fact that the final-state 
jets in pp come predominantly from gluons, as in 
resolved photoproduction. From these compar- 
isons and the ones shown before, it can be con- 
cluded that the pattern of QCD radiation within 
a quark jet is to a large extent independent of the 
hard scattering process. 

4. Quark- and gluon-jet properties 

Several analyses were done at HERA in which 
quark and gluon jets were isolated and the dy- 
namics of the underlying subprocesses studied. In 
some of these analyses, quark and gluon jets were 
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Figure 10. Mean integrated jet shape in inclusive- 
jet and dijet in photoproduction and inclusive-jet 
in NC DIS. 



identified by using samples of dijet events in pho- 
toproduction which contain charm quarks in the 
hnal state. In these analyses, one of the jets is 
tagged as a charm and the substructure of the 
other ("untagged") jet was studied. This method 
provides an unbiased sample of quark jets. In 
direct photoproduction, the untagged jet is also 
a charm quark, but in resolved processes, since 
there are several contributing processes, such as 
gluon-gluon fusion or charm-excitation processes, 
the untagged jet can be a quark or a gluon. 
The ZEUS Collaboration identified quark and 
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Figure 11. Mean integrated jet shape in inclusive- 
jet in photoproduction and NC DIS. 
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jet in NC and CC DIS. 



gluon jets by selecting a sample of dijet events in 
photoproduction and tagging the subsample with 
charm production via the identification of a D* 
meson associated with one of the jets [21]. The 
substructure of the untagged jet was measured 
and found to be consistent with the predictions 
for quark-initiated jets (see Fig. 114b). The sub- 
structure of gluon jets was extracted by means 



of O. 



dijet 



fa ■ O, 



quark 



f 9 -o g 



where the 



substructure of dijets is the measured observable, 
the substructure of quarks was approximated us- 
ing the measured substructure of untagged charm 
jets and the fraction of quarks and gluons was 
estimated using the MC predictions, which de- 
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Figure 15. Mean integrated jet shape in charm 
photoproduction. 



scribe the data well. Figure [Tib shows the ex- 
tracted substructure of gluon jets compared with 
the measured substructure of untagged jets. The 
gluon jets are observed to be broader than the 
quark jets. The model predictions describe the 
measurements well. 
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Figure 14. Mean integrated jet shape in charm 
photoproduction. 



The HI Collaboration made a similar analy- 
sis, but tagging the charm jets by associating a 
muon to one of the jets [IT] . The substructure 



of the other jet was studied for x 



obs 



0.75 (see 



Fig. [ISJ). The QCD-based MC predictions de- 
scribe the data well for x° bs > 0.75. However, 
some differences are observed for x° hs < 0.75. 
The data suggest a smaller fraction of gluon jets, 
which can come only from the gluons in the 
charm-excitation subprocess, at low a;° bs than the 
predictions. 

The differences between quark and gluon jets 
were also studied by exploiting the different type 



of parton content in the final state for one-jet 
and two-jet samples in NC DIS in the laboratory 
frame [22]. The one-jet sample is expected to 
be dominated by quark jets, whereas the two-jet 
sample is expected to contain a larger fraction of 
gluon jets. Figure ITTJI shows the measurement of 
the mean integrated jet shape as a function of r 
for > 14 GeV and for 14 < E^ < 17 GeV for 
the one-jet and the two-jet samples. For the two- 
jet sample, the substructure of the lowest- E 1 ^" jet 
was studied for those events in which the two jets 
are closer than two units in the rj~ <j) plane. These 
jets are observed to be broader than the jets in 
the one-jet sample. The data are compared to 
NLO calculations [13123] . The calculations give a 
good description of the data and show that the 
measurements are consistent with a higher gluon 
content in the two-jet sample. 
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Figure 16. Mean integrated jet shape in one-jet 
and two-jet in NC DIS. 



Since the QCD prediction that gluon jets are 
broader than quark jets has been amply proven 
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by data, it is possible to exploit this fact for iden- 
tifying the different type of jets and study the 
dynamics of the subprocesses. Quark and gluon 
jets were identified [TU] on a statistical basis by 
selecting "narrow" and "broad" jets by requiring 
that the measured jet shape is smaller than 0.6 
or larger than 0.8, respectively. Using this selec- 
tion, the inclusive-jet cross section in photopro- 
duction was measured function of if et (see 
Fig. \T7\ . The measured cross section for broad 
jets displays a very different shape than that of 
the narrow jet sample. The if et distribution for 
the narrow jet sample peaks at about 0.5, whereas 
the broad jet sample distribution increases as ry jet 
increases. The predictions from the QCD-based 
MC models, with the same jet shape selection, 
give a good description of the shape of the data. 
The MC predicts that the broad jet sample is 
dominated by qg — > qg via gluon exchange in re- 
solved processes, whereas the narrow jet sample 
is dominated by boson-gluon fusion events in di- 
rect processes. Therefore, the narrow jet sample 
is indeed dominated by quark-initiated jets in the 
final state and the broad jet sample has a higher 
content of gluon-initiated jets. The comparison 
with the predictions for samples of pure quark or 
gluon jets supports this conclusion. 

ZEUS, 

• ZEUS 98-00 7p (broad jets) 
o ZEUS 98-00 yp (narrow jets) 
— PYTHIA 

3 350 
300 




and gluon exchange. Samples selected accord- 
ing to the internal structure of the jets give an 
unbiased handle to explore the dynamics of the 
subprocesses. Figure \TEk shows the dijet cross 
section as a function of cos((9*) for samples of 
two broad or two narrow jets. The cross sec- 
tion for two broad jets rises more steeply than the 
sample of two narrow jets. The MC predictions 
give a reasonable description of the data. The 
different slope observed in the data can be un- 
derstood in terms of the dominant subprocesses 
in the two samples: the two broad jet sample 
is dominated by subprocesses mediated by gluon 
exchange whereas the two narrow jet sample is 
dominated by subprocesses mediated by quark 
exchange. The cos(#*) distribution for a sample 
of one narrow and one broad jet in the same event 
shows a clear asymmetry (see Fig.fTSb). The MC 
reproduces the data well. The asymmetry ob- 
served can be understood in terms of the dom- 
inant subprocess: qg — ► qg, the positive side is 
dominated by i-channel gluon exchange, whereas 
the negative side is dominated by u-channel quark 
exchange. 
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Figure 18. Dijet cross sections in photoproduc- 
tion. 



Figure 17. Inclusive-jet cross sections in photo- 
production. 



The underlying parton dynamics is reflected in 
the distribution of the scattering angle in the dijet 
centre-of-mass system, 9*. The cos(#*) distribu- 
tion is sensitive to the spin of the exchanged par- 
ticle: a different behaviour is expected for quark 



5. The pattern of parton radiation 

The pattern of parton radiation within a jet 
was studied using sub jet distributions. The pat- 
tern of parton radiation from a primary parton 
is dictated in QCD by the splitting functions, so 
the measurements provide a direct test of these 
functions and their scale dependence. The subjet 
distributions can also be used to study colour- 
coherence effects between the initial and final 
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states. QCD predicts that soft gluon radiation 
tends to be emitted towards the proton direction. 
The measurements of sub jet distributions were 
done in NC DIS [24 as functions of the fraction 



of subjet transverse energy, 
ference between the 77 and 



E: 



sbj /jB jet 



the dif- 



of the subjet with 



respect to the jet, rj sh i~if ct and (/> sbj — <^ et |, and 
a sb j, the angle, as viewed from the jet centre, 
between the subjet with higher- Et and the pro- 
ton beam line in the 77 — <f> plane. The cross sec- 
tions were measured for those jets which contain 
exactly two subjets for y cut = 0.05. Figure fT9l 
shows the normalised subjet cross sections. The 
measurements show that the two subjets tend to 
have similar transverse energy. The distribution 
in rf°i—rf et has a two-peak structure. The cross 
section as a function of \(ff h ^ — <^ et | shows a sup- 
pression around A<f) = 0, which comes from the 
fact that the two subjets cannot be resolved when 
they are too close together; this is also observed 
in the dip between the two peaks in f7 sb j —r] iet . 
The a sh) distribution shows that the highest-i/p 
subjet tends to be in the rear direction. This is 
consistent with the asymmetric two peak struc- 
ture observed in the r] sh i — T] iet distribution and 
related to colour-coherence effects. NLO QCD 
predictions [13j are compared to the data and give 
an adequate description. 
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Figure 19. Normalised subjet cross sections in 
inclusive-jet NC DIS. 



To study in more detail the colour-coherence 



was measured for those jets with subjets of dif- 
ferent transverse energy fraction, separately for 
the low- and high- Et subjets. The measurements 
show that the high-ET subjet tends to be in the 
rear whereas the \ow-Et subjet tends to be emit- 
ted in the forward direction (see Fig. [20]) . The 
NLO predictions describe the data well. This 
behaviour can be attributed to colour-coherence 
effects between the initial and final states, and 
indicates that soft gluon radiation is emitted pre- 
dominantly towards the proton direction. 
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Figure 20. Normalised subjet cross sections in 
inclusive-jet NC DIS. 



The pattern of parton radiation was further 
studied by comparing the predictions for quark- 
and gluon-induced splittings separately with the 
data. The NLO calculations predict that the rate 
is dominated by quark-induced processes. The 
predictions of these two types of processes have a 
different shape: the quark-induced processes have 
subjets with more similar transverse energy than 
the subjets from the gluon-induced processes, and 
the two subjets coming from a qg pair are closer 
together than in the case of qq pairs. The com- 
parison with the measurements shows that the 
data are better described by the calculations for 
jets arising from the splitting of a quark into a qg 
pair (see Fig. [2Tj). 

6. Summary 

Jet substructure has been extensively studied 
at HERA in terms of jet shapes and subjet mul- 
tiplicities and distributions. The measurements 
performed allowed stringent tests of perturbative 
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Figure 21. Normalised sub jet cross sections in 
inclusive-jet NC DIS. 



QCD directly beyond leading order, comparison 
of the properties of quark and gluon jets, the com- 
parison of the pattern of QCD radiation in differ- 
ent hard scattering processes, determinations of 
the strong coupling, the study of the dynamics of 
the underlying subprocesses and the study of the 
pattern of parton radiation. Therefore, jet sub- 
structure provides a powerful tool to make strin- 
gent tests of perturbative QCD. 
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